data_summary<-function(data,varname,groupnames){
  require(plyr)
  summary_func<-function(x,col){
    c(mean=mean(x[[col]],na.rm = TRUE),sd=sd(x[[col]],na.rm = TRUE))
  }
  data_sum<-ddply(data,groupnames,.fun=summary_func,varname)
  data_sum<-rename(data_sum,c("mean"=varname))
  return(data_sum)}

##与fNIRS分析相同，但我有一点疑惑：
##双因素重复检验方差分析中，FR&LR，PRE&POST，均不相互独立，两个因素都是组内因素，属于两因素被试内设计。
fit <- aov(fes_data$bell_error~fes_data$task*fes_data$eva+Error(fes_data$subject/fes_data$eva))
fit <- aov(fes_data$hand_endpoint~fes_data$task*fes_data$eva+Error(fes_data$subject/fes_data$eva))
fit <- aov(fes_data$duration_time~fes_data$task*fes_data$eva+Error(fes_data$subject/fes_data$eva))
fit <- aov(fes_data$shoulder_flex~fes_data$task*fes_data$eva+Error(fes_data$subject/fes_data$eva))
fit <- aov(fes_data$elbow_flex~fes_data$task*fes_data$eva+Error(fes_data$subject/fes_data$eva))
fit <- aov(fes_data$sim2_vector~fes_data$task*fes_data$eva+Error(fes_data$subject/fes_data$eva))
fit <- aov(fes_data$sim2_timeprofile~fes_data$task*fes_data$eva+Error(fes_data$subject/fes_data$eva))
fit <- aov(fes_data$sim2_combine~fes_data$task*fes_data$eva+Error(fes_data$subject/fes_data$eva))
summary(fit)


##配对样本t检验
pre <- subset(fes_data, eva == "pre")
post <- subset(fes_data, eva == "post")
pre_lr <- subset(pre, task == "lr")
pre_fr <- subset(pre, task == "fr")
post_lr <- subset(post, task == "lr")
post_fr <- subset(post, task == "fr")
t.test(pre_lr$pkvel,pre_fr$pkvel,paired = TRUE)
t.test(pre_lr$pkvel,post_lr$pkvel,paired = TRUE)

